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Applying a likelihood analysis to the next-to-mininial supergravity-motivated model, we identify 
parameter space regions preferred by present experimental limits from collider, astrophysical, and 
low energy measurements. We then show that favored regions are amenable to detection by a 
combination of the CERN Large Hadron Collider and an upgraded Cryogenic Dark Matter Search, 
provided that the more than three sigma discrepancy in the difference of the experimental and the 
standard theoretical values of the anomalous magnetic moment of the muon prevails in the future. 
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I. INTRODUCTION 

Supersymmetry is very successful in solving outstand- 
ing problems of the standard model (SM) of elementary 
particles. The theory naturally explains the dynamics 
of electroweak symmetry breaking while preserving the 
hierarchy of fundamental energy scales, it incorporates 
dark matter and the asymmetry between baryons and 
antibaryons, it reconciles the unification of gauge forces 
and accommodates gravity, and more Thus, it is 

important to examine the prospects of the CERN Large 
Hadron Collider (LHC) finding supersymmetry. 

One of the main motivations for supersymmetry is that 
it can naturally bridge the hierarchy between the weak 
and Planck scales. Unfortunately, the presence of the su- 
perpotential fi term in the minimal supersymmetric ex- 
tension of the standard model (MSSM) undermines this 
very aim Experimental data have also squeezed the 
MSSM into fine-tuned regions, creating the supersym- 
metric little hierarchy problem [48]. Extensions of the 
MSSM by gauge singlet superfields not only resolve theu 
problem, but can also reduce the little hierarchy @, H, |j| ■ 
In the next-to-minimal MSSM (NMSSM), the n term is 
dynamically generated and no dimensionful parameters 
are introduced in the superpotential (other than the vac- 
uum expectation values that are all naturally weak scale), 
makingjthe NMSSM a truly natural model MM, S H, B 

E ini, 111, E Q E H [13, H H H i^^ 

Over the last two decades, due to its simplicity and 
elegance, the constrained MSSM (CMSSM) and the min- 
imal supergravity-motivated (mSuGra) model became 
a standard in supersymmetry phenomenology. Guided 
by this, within the NMSSM, we impose the univer- 
sality of sparticle masses, gaugino masses, and trilin- 
ear couplings at the grand unification theory (GUT) 
scale, thereby defining the next-to-minimal supergravity- 
motivated (NmSuGra) model. This approach ensures 
that all dimensionful parameters of the NMSSM scalar 
potential also naturally arise from supersymmetry break- 
ing in a minimal fashion. NmSuGra also reduces the 
electroweak and dark matter fine-tunings of mSuGra. 

Using a simple likelihood analysis, first we identify 
the parameter regions of the NmSuGra model that are 
preferred by the present experimental limits from col- 



lider, astrophysical, and various low energy measure- 
ments. We combine theoretical exclusions with lim- 
its from the CERN Large Electron-Positron (LEP) col- 
lider, the Fermilab Tevatron, NASA's Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) satellite (and other 
related astrophysical measurements), the Soudan Cryo- 
genic Dark Matter Search (CDMS), the Brookhaven 
Muon g— 2 Experiment, and various b-physics measure- 
ments including the rare decay branching fractions b — > 
sj and Bs — > 

Next we show that the favored parameter space can be 
detected by a combination of the LHC and an upgraded 
CDMS, provided that the 3 a discrepancy between the 
experimental and the standard theoretical values of the 
anomalous magnetic moment of the muon (Aa^) per- 
sists in the future. According to the latest calculations, 
the theoretical uncertainty of Aa^ is now under control, 
which makes it a powerful discriminator at the 99 % con- 
fidence level (CL) [H, [24|] [49,] . If, due to new experimen- 
tal input, the Aa^ discrepancy decreases, our results in- 
dicate that allowed regions of the NmSuGra model may 
not be reachable by either the LHC or a 1-ton equivalent 
of CDMS. The detection of these regions will provide a 
new experimental challenge. 



II. THE NMSUGRA MODEL 



In this work, we adopt the superpotential 
W^Wy + XSHu ■ Hd + 



(1) 



where Wy is the MSSM Yukawa superpotential [25|, S 
{Hu,d) is a standard gauge singlet {SU{2)l doublet) chi- 
ral superfield, A and k are dimensionless couplings, and 
Hu ■ Hd = eafjH^H^ with the fully antisymmetric tensor 
normalized as en — 1. The corresponding soft supersym- 
metry breaking terms are 



soft 



-MSSM 

{XA),SHu ■Hd+''^ 



ml\S\^ + 
-S^ + h.c), 



(2) 



-tsoft 
-MSS 

not the term. 



where Cf^/gg]^.[ contains the mass and Yukawa terms but 
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The superpotential ((T]) possesses a global Z3 symmetry 
which is broken during the electroweak phase transition 
in the early universe. The resulting domain walls should 
disappear before nucleosynthesis; however breaking 
(via singlet tadpoles) leads to a vacuum expectation value 
(vev) for the singlet that is much larger than the elec- 
troweak scale. Thus the requirement of the fast dis- 
appearance of the domain walls appears to destabilize 
the hierarchy of vevs in the NMSSM. Fortunately, in 
Ref. [H, [l^] is was shown that, by imposing a R- 
symmetry, both the domain wall and the stability prob- 
lems can be eliminated. Following 26], wc assume that 
tadpoles are induced, but they arc small and their effect 
on the phenomenology is negligible. 

We also assume that the soft masses of the gauginos 
unify to Mi 12-, those of the sfermions and Higgses to Afo, 
and all the trilinear couplings (including and A\) to 
Aq at the GUT scale. This leaves nine free parameters: 
Mo, Ml/2, ^, (^«), (^fd), (^> (the Higgs and singlet 
vevs), 7715, A and k. The three minimization equations 
for the Higgs potential [H and {HuY ^ {Hd)"^ = v"^ (here 
V = g2 is the standard Higgs vev) eliminate four 
parameters. With the introduction of ^ = ^{S), and 
tan (3 = {H^) / {Hj) , our remaining free parameters are 



Mo, Ml/2, ^0, tan/3. A, sign(/i). 



(3) 



Constrained versions of the NMSSM have been studied 
in the recent literature. The most constrained version is 
the cNMSSM Hi] with ms = Mq. In other cases the 
Ak = A\ relation is relaxed [1^, and/or k is taken as 
a free parameter [s^, |3l|, or the soft Higgs masses are 
allowed to deviate from Mo [12] giving less constrained 
models. In the spirit of the CMSSM/mSuGra, we adhere 
to universality and use only A to parametrize the singlet 
sector. This way, we keep all the attractive features of the 
CMSSM/mSuGra while the minimal extension alleviates 
problems rooted in the MSSM. 

Our goal is to show that the experimentally favored re- 
gion of the NmSuGra model can be discovered by nascent 
experiments in the near future. To this end, we use the 
publicly available computer code NMSPEC [3^ to calcu- 
late the spectrum of the superpartner masses and their 
physical couplings from the model parameters given in 
Eq. (O. Then, we use NMSSMTools 1.2.5, an extensively 
modified version of ISATools and that of DarkSUSY 4.1 
[34j to calculate the relic density of neutralinos (f^/i^), 
the spin-independent neutralino-proton elastic scattering 
cross section {a si), the NmSuGra contribution to the 
anomalous magnetic moment of the muon (^QNmSuGra^^ 
and various b-physics related quantities. 

For each set of the model parameters, we quantify the 
experimental preference in terms of 



i:(- 



experiment 



m: 



NmSuGra 



1/2 



(4) 



where is the central value of a physical quantity mea- 
sured by an experiment or calculated in the NmSuGra 
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FIG. 1; Square-root of the exponent of the likelihood function 
vs. four of the NmSuGra input parameters for a random sam- 
ple of models. The combination of the experimental quantities 
included in favor low values of Mo, M1/2 and \Ao\. 



model, and ai is the combined experimental and (where 
available) theoretical uncertainty. The sum includes the 
experimental upper limits [50|for 

(1) Qh^ = 0.1143 ± 0.0034 [35], 

(2) Br{Bs ->/!+//-) = 5.8 x 10"^ (95 % CL) [36|, and 

(3) as I by CDMS 37], 

the LEP lower limits of the lightest scalar Higgs and 
chargino masses (which can be approximately stated as) 

(4) TUh > 114.4 GeV for tan/3 - 10 

TUh > 91 GeV for tan/3 - 10, 

(5) m^^ > 104 GeV, 
and the central values of 

(6) Aa^ = 29.5 ± 8.8 x 10"^° [H, 

(7) Brib ^ 57) = 3.55 ± 0.26 x lO"'* [s^. 

The related uncertainties are given above at 68 % CL, un- 
less stated otherwise. Theoretical uncertainties are calcu- 
lated using NMSSMTools for the b-physics related quan- 
tities. Among the standard input parameters, mi,(mi,) = 
4.214 GeV and mf'^^ = 171.4 GeV are used. 

A glance at experimental likelihood reveals a signifi- 
cant statistical preference for relatively narrow intervals 
of Mo, Ml/2 and ^o- This is illustrated in FIG. [TJ which 



shows vx^ as the function of Mo, M1/2, ^0 and tan/3 
for sign(/i) > for a randomly selected set of models 
[sH . From FIG. [TJ we conclude that it is enough to ex- 
amine the low Mo, Ml/2 and \Aa\ region, since the rest is 
disfavored by the above combination of the experimen- 
tal data at (or more than) 99 % CL. Similarly to the 
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CMSSM/mSuGra, at high values of Mq, M1/2 and |^o|, 
^ is dominated by Aa^. 

Based on the above, we hmit our study to the following 
ranges of the continuous parameters [5^ : 

< A/o < 4 TeV, < M1/2 < 2 TeV, 
0<|Ao|<5TeV, 1 < tan/3 < 60, O.OK A < 0.7. (5) 

The upper limit on A arises from the requirement that its 
running remains perturbative as it evolves up to the GUT 
scale. Restricting Mq, M1/2 and Aq to such values also 
app ears consistent with the electroweak precision data 
[39| and greatly reduces fine tuning. 



III. DETECTABILITY OF NMSUGRA 

Having defined the NmSuGra model and constrain- 
ing the range of its parameters, we set out to show 
that this parameter region will be detectable by the 
LHC and an upgraded CDMS detector. To this end, 
we randomly select about 20 million models in the range 
defined in Eq. (O, and for each model point evaluate 
Eq. (HI). Two million theoretically allowed representa- 
tive model points are projected in FI G, p to the fi/i^ vs. 
{Nl^ + iVf2)/(l - ^11 - ^12) plane This plane of 

f2/i^ and the neutralino mixing matrix elements (J^ij) is 
a good indicator of the gaugino, higgsino and singlino 
admixture of the lightest neutralino. From FIG. [Hit is 
evident that the WMAP upper limit (green horizontal 
line) allows models with mostly bino- (red) and higgsino- 
like (magenta) lightest neutralino, while the fraction of 
allowed models with singlino-like (blue) dark matter is 
negligible. As such the NmSuGra model is very similar 
to the CMSSM/mSuGra. 

The similarity is even more evident when we exam- 
ine the grouping of the WMAP allowed model points. 
By checking mass relations and couplings, we can eas- 
ily establish that the various branches (denoted by 
1, 2, 3 and 4) belong to models with distinct neu- 
tralino (co)annihilation mechanisms well known from 
the CMSSM/mSuGra [4l|. Branch 4 contains only 
models with dominant neutralino-stop coannihilation, 
while branch 3 corresponds neutralino-stau coannihila- 
tion. Branch 2 represents the Higgs resonance corridors, 
and branch 1 is the equivalent of the CMSSM/mSuGra 
focus point region. 

To gauge the detectability of the NmSuGra model, first 
we identify model points that could have been seen at 
LEP. We require the lighter chargino to be lighter than 
103.5 GeV or the lightest scalar Higgs to be lighter than 
114.4 GeV. We relax the latter to mu when the mass of 
the lightest pseudo-scalar Higgs [mn^^ ) approaches that 
of the lightest scalar {mH^^) H^]. We do not apply the 
114.4 GeV LEP Higgs limit when we encounter either 
a mostly singlet lightest Higgs or when mn,^ < 'niHp-^ 
[3^ . [43] . Focusing our attention to the lower right corner 
of FIG. m in the top left frame of FIG. Owe show the same 
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FIG. 2: Relic abundance of the lightest neutralino as the 
function of its gaugino admixture. Right of the vertical line 
the neutralino is mostly gaugino. The horizontal lines shows 
the WMAP upper limit (95 % CL). 



model points in different coloring. Models detectable by 
LEP are colored green, that is the green points repre- 
sent the reach of LEP for the NmSuGra model. The 
red colored models pass the above LEP constraints, i.e. 
are allowed by LEP. Just as in the CMSSM/mSuGra the 
neutralino-stop coannihilation region is mostly covered 
by LEP. 

To estimate the LHC reach, we rely on the similarity 
between the mSuGra and NmSuGra models. Accord- 
ing to Ref. [31 the reach of the LHC for mSuGra can be 
well approximated by the combined reach for gluinos and 
squarks. Based on this, if either the gluino mass is be- 
low 1.75 TeV, or the geometric mean of the stop masses 
is below 2 TeV for a given model point, we consider it 
discoverable at the LHC. While this is an approximate 
statement that has to be supported by a detailed study 
of LHC event generation in the given models, in the light 
of Ref. [44] it is a conservative estimate of the LHC reach. 

The top right frame of FIG. [3] shows the model points 
that can be reached by LEP and the LHC using the above 
criteria. As in the CMSSM/mSuGra, most of the slep- 
ton coannihilation and the bulk of the Higgs resonance 
branches are covered by the LHC. A good part of the 
focus point is also within reach of the LHC, with the 
exception of models with high A/q and/or Mi/2- 

The bottom left frame of FIG. [3] shows the reach of 
a one ton equivalent of CDMS (CDMSIT). As expected 
from the CMSSM/mSuGra, the rest of the focus point 
and most of the remaining Higgs resonances are in the 
reach of CDMSIT. The small number of models that re- 
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(Nj]+Njl)/(1-Nj]-N^l) 

FIG. 3: As FIG. [2] except green models can be reached by the 
combination of LEP, the LHC and CDMSIT, while red ones 
cannot. The last frame dismisses experimentally inaccessible 
points which have ^/x^ > 3. 



main inaccessible are all located in regions that have rel- 
atively low Afo and high Mi with dominant neutralino 
annihilation via s-channel Higgs resonances. The Nm- 
SuGra contribution to Aa^ in these model points is out- 
side the preferred 99 % CL region as shown by the last 
frame. 

While we focused our discussion and plotting on the 
higgsino and gaugino regions, we carefully checked that 
LEP, the LHC and the upgraded CDMS experiments 
combined can also detect models with WMAP allowed 
singlino-like dark matter. The reason for this is the fol- 
lowing. The singlino-like lightest neutralinos typically 
take part in more than one (co) annihilation mechanisms 
to satisfy the WMAP limit. While being on a Higgs res- 
onance, they also coannihilate with sfermions. In this 



case, either being in the sfermion coannihilation region 
ensures the LHC detectability or the Higgs resonance is 
strong enough to enhance recoil detection. 

Furthermore, the universal mass relations ensure that 
the lightest scalar Higgs never decays to a pair of the 
pseudo-scalar Higgs bosons in the phenomenologically al- 
lowed region of NmSuGra, just as in the cNMSSM ^29j. 
Similarly, the lightest singlinos, allowed by the present 
experimental constraints are always heavier than about 
100 GeV. Assuming that the NmSuGra contribution to 
the anomalous magnetic moment of the muon is larger 
than a minute 3.1 x lO""'^" constrains slepton and chargino 
masses below 3 and 2.5 TeV, respectively. Since uni- 
versality restricts the mass hierarchy within NmSuGra, 
the resulting mass spectrum is typically mSuGra-like. 
Thus, the cascade decays and their signatures at LHC 
are not expected to be significantly deviate from that 
of the mSuGra case. The most typical NmSuGra decay 
cascade at the LHC would be giuino — )■ squark, quark — > 
chargino/neutralino, W/Z — + neutralino, SM particles. 

IV. CONCLUSIONS 

Analyzing the next-to- minimal supergravity motivated 
(NmSuGra) model, we found that the LHC and an up- 
graded CDMS experiment will be able to discover the ex- 
perimentally favored region of this model, provided that 
the present deviation between the experimental and stan- 
dard theoretical values of the muon anomalous magnetic 
moment prevails. If, due to future experimental input, 
the (? — 2 constraint weakens, then certain parameter re- 
gions of the NmSuGra model will have to be detected in 
alternative ways. 
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